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The telemetry performance of an arrayed receiver system, including radio losses, is 
often given by a family of curves giving bit error rate vs bit SNR, with tracking loop SNR 
at one receiver held constant along each curve. This study shows how to process this 
information into a more compact, useful format in which the minimal total signal power 
and optimal carrier suppression, for a given fixed bit error rate, are plotted w data rate. 
Examples for baseband-only combining are given. When appropriate dimensionless vari- 
ables are used for plotting, receiver arrays with different numbers of antennas and 
different threshold tracking loop bandwidths look much alike, and a universal curve for 
optimal carrier suppression emerges. 


I. Introduction 

In an internal memorandum, J. W. Layland wrote: 

The task of comparing the performance of the 
Arrayed Deep Space Network with and without 
combined carrier references (CCR) seems to be 
needlessly complex. While the analysis machinery 
which has been built for us is capable of precisely 
describing the link performance under a plethora 
of conditions, interpretation of those results in 
terms useful for deciding whether to deploy the 
CCR or not is difficult, and potentially even 
ambiguous. There are too many “free” variables 
which seem to be available — even though physical 
constraints will interrelate many of them. Our task 
is somewhat analogous to driving a car with a 
separate steering wheel for each front wheel. 

For example, the telemetry performance of a DSN receiver 
array can be represented by a family of “waterfall” curves 
giving bit error rate as a function of array total Ej^INq , the 


sum of the bit signal-to-noise ratios at all the array elements, 
with p^j, the tracking loop SNR at the strongest array ele- 
ment, held constant (Refs. 1 and 2 and Fig. 1). We need such a 
plot (a family of six curves in this instance) for each channel 
code, each DSN antenna array, each signal combining method, 
and each choice of threshold tracking loop bandwidth w^q = 
2Z>^q . (It is assumed that all receivers in the array use the same 

*^to) 

This presentation is more convenient for the person gener- 
ating the plots than for the person who has to interpret them. 
To reduce the volume of results, one can settle on a fixed 
telemetry performance, Pg = 0.005, for example, and use each 
waterfall curve set to generate one “radio loss” curve: 

V^o‘(V^o)oo vsp^, , 

where £’j/Nq is the array total bit SNR needed to attain Pg 
for the given p^j , and {EJNq)„^ is the SNR needed to attain 
Pg for p^j = + °°, i.e., perfect phase tracking (Ref. 1). 
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In the memorandum quoted above, Layland suggested 
another presentation format, in which one plots total power 
threshold as a function of data rate, again for a fixed P^. 
Further, for each. data rate, one uses the optimal modulation 
angle 6, i.e, the 6 that requires the least total power. This 
eliminates 0 as a free parameter and tells how to design the 
signal. 

In the present paper, this idea is applied to a set of results 
for baseband-only combining. As Layland suggests, we hold to 
a bit error rate of 0.005. It turns out that a judicious choice of 
dimensionless variables allows us to use a “universal” curve for 
the optimal 6 (or carrier suppression cos^fl) as a function of 
data rate. Further, when regarded from this point of view, all 
the arrays look the same within 0.5 dB. In other words, we 
have a nearly universal curve for total power-to-noise ratio vs 
data rate. The differences among the arrays can then be seen 
by plotting EJNq vs data rate. 

II. A Sample of Results 

The source data come from internal memoranda of 
D. Hansen and D. Divsalar, and from unpublished results of 
L. Deutsch. Let us list the conditions: 

(1) (7, 1/2) convolutional code, maximum likelihood 
decoding. 

(2) One-way radio losses only. 

(3) High-rate bit error probability model. 

(4) Baseband-only combining with weights optimal for per- 
fect phase tracking. 

(5) = 10 Hz or 30 Hz, the same for all receivers in the 
array. 

(6) Antenna array combinations (diameters in meters): 64, 
64-34, 64-34-34, 64-34-34-34. We assume that gain-to- 
noise ratio is proportional to area. 

For each array and there is a set of six waterfall curves 
giving /*£. vs EjNg forp^j = 10, 11, 12, 13.5, 15, and 40 dB. 
Figure 1 is the plot for the four-element array and = 
10 Hz. 

Setting a telemetry performance threshold of Pg = 0.005, 
we compute optimal modulation angle and minimal total 
power-to-noise ratio. Results are plotted in Figs. 2-4. 

The dimensionless variables used for plotting are listed 
below. In the plots, the term X (dB) always means 10 logjQ X. 

(1) = normalized data rate, where = data rate. 

(2) cos^0 = carrier suppression. 


(3) Pj’KNqWj^q) = normalized array total power-to-noise 
ratio, where PjiNg is the sum of the P-j-lNg values of 
the receivers in the array {Pj, = carrier power + data 
power, Nq = 1 -sided noise spectral density). 

(4) E^lNg = array total bit SNR, the sum of the Ej^INg 
values of the receivers in the array (£’j = data energy 
per bit). 

Figure 2 is a “universal” curve giving optimal carrier sup- 
pression vs normalized data rate. The modulation angle is 
truncated at 80 deg. Actually, the several arrays yield distinct 
curves, which all agree within about 1 dB. Since the optimiza- 
tion problem is smooth, a deviation from optimal carrier 
suppression causes only a second-order increase in the total 
power needed to maintain the required performance. This is 
why one curve will do for all the array- combinations. 

Figure 3 shows the array total Pj^/Ng needed to achieve P^ 
= 0.005. Total power has been minimized by the proper choice 
of carrier suppression (Fig. 2). Although this curve is again 
billed as universal, it is actually about 0.5 dB thick in the 
vertical direction. Thus it is intended only as a long-range 
snapshot. We shall soon tell how to compute values with 
greater accuracy. 

For each array, Fig. 4 shows the array total £),/Aq needed 
to achieve P^ = 0.005. Since Ef^lNg varies much less than 
PtII^O’ "O"' easily distinguish the different arrays. 

Still, the curves all fit within a 0.5-dB band. Having these 
EbINg curves, one can compute Pj’INg more accurately than 
before by the formula 


1 - COS^ 0 -^0 ^LO 

where cos^0 is read from Fig. 2 and EJN^ from Fig. 4. 


III. Remarks and Cautions 

These results, though incomplete and inaccurate, are ade- 
quate for a pilot study. The E^/Ng curves are probably good 
within ±0.1 dB. The author’s source data consist of meager 
sets of numbers read from photocopied computer-generated 
graphs. Although both the Hansen-Divsalar and the Deutsch 
graphs were processed, we show only the Hansen-Divsalar 
results because Deutsch’s basic Viterbi decoding error curve 
(used in the high-rate model) is slightly different, and his 
graphs are harder to read accurately. Thus, for each array there 
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is only one . The Deutsch results indicate that the Ef^ /Nq 
curves for vv^q = 10 Hz lie perhaps 0.1 dB above the corre- 
sponding curves for = 30 Hz. This, of course, does not 
mean that narrowing the loop bandwidth makes telemetry 
worse; recall that data rate is normalized by Wj^q . Thus for a 
given array, if we were to plot the 10 Hz and 30 Hz E^INq 
values against unnormalized data rate, the 10-Hz curve would 
lie to the left of the 30-Hz curve, and therefore below it. 

Notice, too, that the EjN^^ curve for any multiple-antenna 
array lies above the curve for one antenna. This is because we 
are, in effect, comparing the array to a single large fictitious 
antenna with the same total area. The array is worse because 
the tracking loop of the fictitious receiver gets more carrier 
power than any of the loops of the array receivers. To com- 
pare a 64-m station to the 4-element array on the basis of 
£’j,/jVq at the 64-m station, one can simply move the 4-element 
curve down by 1 + 3(34)^/64^ = 2.66 dB. 

A particular point on any of these curves is valid only if 
Pj^j , the loop SNR at the strongest array element, is at least 
10 dB. The left-hand endpoint.of each EJNq curve shows the 
value beyond which this condition fails. For the array 


in question, the “universal” carrier suppression and total 
power curves must not be used to the left of this value. 

IV. Conclusions 

Layland’s suggestion for presenting array telemetry per- 
formance appears to work well when applied to a set of results 
for baseband-only combining. Each set of waterfall curves is 
distilled into one curve of EJNq vs normalized data rate 
Pe constant and modulation angle opti- 
mized. By plotting array total Ef^lN^ and using as a 
normalizing factor, we absorb the gross differences among the 
different antenna arrays and receiver bandwidths, thus leaving 
a set of curves that differ by at most 0.5 dB. Further, there 
emerge approximate universal curves for optimal carrier sup- 
pression and minimal array total P j,I{NqWj^q) 

Although this work is merely a small pilot study, the 
method may be useful for comparing a variety of situations, 
including various carrier combining methods and channel 
codes. If this is the case; then the method could also be used 
for handbook specification of the arrayed DSN. The carrier 
suppression curve (Fig. 2) may be convenient for mission 
design. 
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Fig. 1. A family of waterfall curves giving bit error rate vs array 
total with = tracking loop SNR at the largest antenna. 

Conditions: (7, 1/2) code, one-way radio losses, baseband-only 
combining, 10-Hz tracking loops, four-element array: 64-34- 
34-34 m (after Hansen) 



Fig. 2. An approximate universal curve for optimal carrier 
suppression vs data rate normalized by threshold loop bandwidth. 
Conditions: (7, 1/2) code, one-way radio losses, one to four DSN 
stations, baseband-only combining, P^ = 0.005. The modulation 
angle 0 is truncated at 80 deg 





Fig. 3. An approximate universal curve giving normalized array 
total power-to-noise ratio vs normalized data rate, for an error rate 
Pf = 0.005. Carrier suppressions from Fig. 2 are used. This curve 
actually consists of several curves differing by at most 0.5 dB. 
Conditions: same as Fig. 2 





Fig. 4. Array total E^/Afg needed to achieve an error rate P^ = 0.005, 
with carrier suppression from Fig. 2. Arrays: (1) One antenna, w^g = 
30 Hz; (2) 64 m - 34 m, 10 Hz; (3) 64-34-34 m, 10 Hz; 
(4) 64-34-34-34 m, 10 Hz. Conditions: same as Fig. 2 
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I. Tracking Loop Formulas 


= two-sided loop noise bandwidth 


We are given a set of waterfall curves, giving vs array 
total EJNq with fixed at 10, 11, 12, 13.5, 15, and 40 dB. 
Since loop SNR p^ is a nonlinear function of carrier margin m, 
a more well-behaved quantity, we first convert the p^ values to 
the equivalent m values. Recall that the DSN definition of m is 


0.862 + p^ 

= limiter performance factor 


(A-8) 


m 


P 

C 


^O^LO 


where P^ = carrier power. 

For the reader’s convenience, here is a set of formulas 
relating p^ to m for the DSN carrier tracking loops. 


This is a rough approximation due to Tausworthe (Ref. 3). For 
more accuracy, one should use a formula of Springett and 
Simon (Ref. 4): 


1 - exp (-p^) 

a2(l -t 0.0975 exp (-0.2146 p^)) 


(A-9) 


= 2-sided IF predetection bandwidth 


Then we have 




2w 


LO 




H 


(A-1) 


2m w 


P, = 


LO 




m 


(r+DF 


(A-10) 


= threshold predetection SNR 

Define 

(Hx) = e~^[I^{x)-yi^{x)] (A-2) 


This gives loop SNR in the form p^ = To obtain 

m from p^, one can iterate the formula m = 
from some initial guess for m. The iteration, which is conver- 
gent because i// is a slowly varying function, can be accelerated 
by Steffensen’s method (Ref. 5). 


Then 


2 "^^PhO ^{2 ^Ho) 


li. Optimal Carrier Suppression 

(A-3) 

Having converted each pj^j to m^ = carrier margin at the 
strongest array element, we multiply by array total 


= threshold limiter suppression factor 


antenna area/area of largest antenna to obtain total carrier 
margin m. We draw a horizontal line across the set of waterfall 

= '^Pho 

(A-4) 

curves at the desired P^ (5 X 10“^ in this study). The 
intersections of this line with the waterfall curves give several 

= predetection SNR 


points [m^g, (£"fc/A^o)<is] > through which we interpolate a 
smooth function 


(A-5) 

(^-11) 


= limiter suppression factor 
r = 2oc/a^ 


by means of cubic splines or a similar method. (The author 
used a local cubic method of Butland and Brodlie, Ref. 6.) 
(A-6) Then we can also write 


w. 


r+ 1 
3 


w 


LO 


(A-7) 


^ = fim) = 
■'''o 


(A- 12) 
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Define the dimensionless variables 


original page ® 

^'^R QUALfTf 


Since pint = sec^ d, we have 


^O^LO 




where P^INq is the array total. Then we have the fundamental 
equation 


tan^ e = -\lg'{m.„) 


(A-17) 


We use Eq. (A-17) to determine 6, reduce it to 80 deg if 
necessary, and compute p and r as before. It is curious that 
optimal modulation angle can be so easily extracted from the 
slope of g. 


p = m + rf{m) 


(A- 13) 


This just says total power = carrier power + data power, and 
data power is set for the desired P^. Now, given r (i.e., data 
rate), we ask for the minimal p (i.e., total power) that gives the 
desired . To get it we just set dpidm = 0: 


0=1+ rf\m) 
fXm) = - Mr 


(A- 14) 


If we solve Eq. (A-14) for m, then we can get p from Eq. 
(A- 13). The optimal carrier suppression ism/p. 

There are two tricks for simplifying this procedure. First, 
there is no need to solve Eq. (A-14) for m in terms of r. 
Indeed, if we use m as independent variable instead of r, then 
Eqs. (A-14) and (A- 13) give r, p, and m/p as functions of m\ 
we then plot p, E/^INq , and m/p vs r. If the modulation angle 6 
(rn/p = cos^ 6) comes out greater than 80 deg, then set 

p = m sec^ (80 deg) 
and solve for r from Eq. (A- 13). 

The second trick comes from using Eq. (A-1 1) directly. By 
Eq.(A-12), 


= mfXm)lf(m) 

which, combined with Eqs. (A-13) and (A-14), gives 

rf(m) 1 
m g'{m^) 


(A- 15) 


III. Adjusting for the Universal Carrier 
Suppression Curve 

Carrying out all of the above for several array configura- 
tions, we may find that the curves of optimal cos^0 vs r almost 
fall on top of each other, within 1 dB, say. We can then 
replace these curves by one universal curve. If, however, we 
combine a universal d with an EJNq from above, then Eq. 
(1), which in our present notation reads 




(A-18) 


gives a wrong answer for p. 


For the user’s sake, it is necessary to apply a small fudge 
factor to E^jINq. As we mentioned, a small deviation of d from 
the optimum hardly changes p, since p has been minimized as 
a function of 6. Given r, we have exact and universal values of 
cos^d, and an exact value of E^Nq. An adjusted value of 
EjjINq is obtained from Eq. (A-18) by pretending that 
adjusted p = exact p: 

(EjNo)^i _ exact 

sin^ 6 . sin^ 0, 


These adjusted values were used for Fig. 4. 


As a check, we can derive a new value of p from 


(E./N^) .. = f(p cos^ 6 . ) 

b O'adj ^new univ' 


^= 1 - 




(A-16) 


The new p might be about 0.01 dB greater than the old p 
derived from inserting either {_EJNq, 6) pair into Eq. (A-18). 
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